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Abstract 

The Tucunduba Dam is about 290 km West of Fortaleza, Ceara 
state. The seismic monitoring of the area, with a local network, has 
began on June 11, 1997, soon after the occurrence of the event with 
magnitude 3.2 on June 09, 1997. The monitoring was done with 
one the analogical station (used for magnitude duration measurement 
and statistical control activity) and seven digital stations. The digi- 
tal stations, with three components each, were operated from June to 
November 1997. In this work, the data collected during digital mon- 
itoring was analyzed to determine hypocenters and focal mechanism. 
To determine hypocenters, the HYP071 program was used, with half- 
space model, whose parameters were 5.95 km/s, for P-wave velocities, 
and 1.69 for the ratio between P and S-wave velocity. The active zone 
was nearly 1 km length, with depth between 4.5 and 5.2 km. With 16 
events recorded in the same six stations, we determined the direction 
of the fault plane (NE-SW). The fault mechanism is strike-slip with 
a small normal component. The dip estimate was 65°SE with FPFIT 
and 80° SE with FOCMEC. Preliminary estimates of maximum hori- 
zontal compressive stress, using P-axis direction, were in accordance 
with Ferreira et al. (1998). The small difference is probably due to 
influence of the sedimentary basin in the regional stress. The active 
area is in accordance with seismicity described by Assumpgao (1998), 
i. e., the combination regional stress, local flexural effect from thick 
sedimentary loads, and a presumably weaker crust from Mesozoic thin- 
ning. 
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1 Introduction 



The Northwest of Ceara is an area of seismic activity known since the 
XIX century (Ferreira & Assumpgao, 1983; Berrocal et al., 1984; Ferreira 
et al, 1998). On June 09, 1997 (district of Serrota, Senador Sa-CE), an 
event with magnitude m& = 3.2 (Intensity III-IV MM) happened, scaring the 
local population. The growth of the seismic activity in Serrota, specially in 
the Tucunduba Reservoir, led the UFRN to study the seismic activity, to 
determine the hypocenter and displacement mechanism with large precision 
and to try to understand the local seismicity. 

For this study, an analogical station and a digital network with seven 
stations were installed (Figure [TJ. The analogical station operated from 
June 11, 1997 to August 14, 1998 and the digital network from June 18 
to November 05, 1997. The analogical station was used in the statistics 
of the earthquakes as well as to obtain the magnitude with the duration, 
using the same parameters presented by Blum & Assumpgao (1990), since 
the analogical station operated in the same conditions. 

2 Geological comments 

The seismic area is part of the Borborema province which is divided in 
five tectonic domains (Brito Neves et al., 2000). The Northwest of Ceara is 
divided by two of these domains; the Medium Coreau Domain (MCO) and 
the Ceara Central Domain (CC). 

The MCO Domain is situated between the margin of the Sao Luis- West 
Africa craton and the Transbrasiliano/Kandi Lineament. The domain con- 
sists of basement (2.5 Ga) high-grade methamorphic rocks and pelitic-carbo- 
nate fold belts. Along the Transbrasiliano lineament occured a series of 
transtensional basin transpressional and post-orogenic plutons occupying 
pull-apart; many of these are covered by the Phanerozoic sedimentary rocks 
of the Parnafba Basin (Gusmao, 1998). 

The CC domain is situated between the Transbrasiliano/Kandi lineament 
and the Senador Pompeu lineament (Figure EJ). This domain consists of 
gneissic basement formed during the Transamazonian collage, with invulsion 
of an Archean nucleus. The basement has 2.2 and 2.1 Ga. The domain 
also contains a series of middle Neoproterozoic supracrustal sequences or 
remanent of folds belts and expressive Brasiliano plutonism. 

The dam is located in the bound of the MCO domain and coverings 
Phanerozoic. Close to the dam, there is a sucession of horst and grabens 
(Torquato & Nogueira Neto, 1996). 



2 



3 Seismic Activity 



The large event recorded in the Tucunduba Dam, with magnitude 3.2 
nib, was on June 9th, 1997 and was felt mainly in Serrota, a small villages 
of Senador Sa, by the margins of the dam (Figure P), where roofs were 
shaken, glasses of stores balanced, school materials were dropped off desks 
and, the inhabitants were afraid (Franga, 1999). Then, on June 11th, 1997, 
the analogical seismograph station, with MEQ-800 smooked paper recorder 
(SN1A), was installed. Seven digital short period portable stations were also 
installed, on June 18, 1997. Except the station SN03, all the stations were 
installed in granitic/gneissic bedrocks of the crystalline basement, allowing 
low noise and clear P and S wave arrivals. The digital stations clocks were 
corrected once each hour with GPS and P-wave readings had accuracy to 
±0.001s. 

The station SN1A operated from June 11, 1997 to August 14, 1998, 
recording a total of 2217 events (Figure OJ), and the largest magnitude was 
m;, = 3.0 on December 12, 1997. 

3.1 Determination of Hypocenters 

Hypocenters were determined with the program HYP071 (Lee & Lahr, 
1975). The half-space velocity model has been producing good results when 
applied to several areas of Brazilian Northeast (Ferreira et ai, 1987; 1995; 
1998; Takeya et ai, 1989; do Nascimento, 1997). This half-space model was 
adopted, since the seismic area is in the precambrian basement composed 
of consolidated rocks and low attenuation, generating clear P and S waves 
arrivals in the seismographs (Figure 0]). 

Various tests were made to find the best P velocity (vp) and (fp/fs) 
ratio. We varied P wave velocity between 5.4 and 6.45 km/s and vp/vs ratio 
between 1.60 and 1.74. The velocity model with the lowest root mean square 
(rms) time residual had vp/vs = 1-69 and vp = 5.95 km/s. The value for vp 
is acceptable, since in different parts of the world, the vp in granitic/gneissic 
upper crustal rocks varies from 5.8 to 6.4 km/s (Christensen, 1982; Chris- 
tensen & Mooney, 1995). For this model, all 160 selected events (recorded 
by at the least four stations) had rms residual 0.03 sec and horizontal errors 
less than 0.5 km. The most of the events have rms < 0.01, erh < 0.1 km 
and erz < 0.1 km, constitutes a set of weel determinated. This is basically 
due the site of deployment (mainly bedrock), the stations near the epicentral 
area anda sample rate of 500 sample per second. 

Figure El shows the epicenters of the events with rms residuais less than 
0.03 sec and horizontal errors less than 0.5 km, recorded during the digital 
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network operation. As can be seen, all events occurred inside the lake, with 
about 1 km length and depths between 4.5 and 5.2 km. We do not know 
the effect of the water in the trigger, since the dam is very old, shallow and 
seismic activity is not known previous 

4 Fault plane 

To estimate the fault plane, we used 16 well-located events recorded at 
the same six stations (Table [TJ Figure EJ. The hypocenters were distributed 
in 0.3 km length, with depth varying from 4.68 to 4.87 km (Figure |7J). 

Starting from these data, we determined, by least squares, the direction 
and the dip of the fault plane. The azimuth was 60° ± 10° and the dip 
88° ± 11°. Projections of the hypocenters used in the calculation, in parallel 
and perpendicular plans to the fault plan, are shown in the Figure The 
vertical and horizontal maximum errors observed for that set of data were 
0.1 km, while the range of the active zone was of 0.3 km and the variation in 
depth was of 0.2 km. In these conditions, the hypocentral data can provide 
only a rough estimate of the fault plane orientation. 

5 Focal mechanisms 

Composite focal mechanisms were determined using clear P wave first mo- 
tions with aid of the FPFIT program (Reasenberg & Oppenheimer, 1985). 
FPFIT program does a grid search to find the solution that minimizes a 
weighted sum of discrepancies in the polarities, considering both the esti- 
mated variance of the data and the theoretical P wave radiation amplitude 
(Reasenberg & Oppenheimer, 1985). The individual focal mechanisms were 
determined using S/P amplitudes ratios in additition to P wave first motions, 
according to the method of Kisslinger (1980) and Kisslinger et al. (1982). For 
individual focal mechanism, we used the FOCMEC program (Snoke et ai, 
1984), which does a grid search of the nodal parameters to find the solution 
the best fits the observed log (amplitude ratio). 

Firstly, with the 16 events of Table Q we used the fault plane obtained 
by least square with reference, we estimated the composite focal mechanism 
(with number of P + S readings > 12 and at least five clear P onset). Added 
the short-period station located in Sobral-CE (SB17), epicentral distance of 
85.5 km, for the two largest events in Tabled (Event 1 and 2). To deter- 
mine the focal mechanism by using FPFIT, we limited the fault plane near 
60° with range 30° and result of the composite focal mechanism (azimuth, 
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dip and the rake), shown in the Table |2hFigure 01 The solution of the focal 
mechanism is approximately a NE-SW fault, with strike-slip dextral move- 
ment with normal component. The fault plane was chosen by taking into 
account the hypocentral distribution (Figure EJ). The composite focal mech- 
anism solution is a strike-slip type, which is predominant in the northeast 
of Brazil (Assumpgao et al., 1985; 1989; Ferreira et al., 1995; 1998). The 
direction of the P axis (278°), obtained by FPFIT, differs from 6° of the 
maximum horizontal compressive stress (SH max ), obtained by Ferreira et al. 
(1998) in NE Ceara region. 

Secondly, seven events of Table [T] were selected (Table EJ), by using as 
selection criterion the largest number of clear polarities of the S wave. The 
results obtained through the FOCMEC are in Table H3 limiting the range 
of orientation and dip of the axes B (null axis), and A and N (correspond 
to Herrmann's X and Y axes), with step 2°. The events are numbered in 
accordance with Table Q and the best solutions of focal mechanisms with its 
P and T axes, for the reliable mechanisms for FOCMEC. The results showed 
a predominance of strike-slip mechanism (1, 11, 14, 15 and 16). Event 6 
has normal mechanism and event 13 a normal with strike-slip component. 
Without considering the dip, several mechanisms have close direction to the 
estimated azimuth with previous methods (60°) (Table 

The average strike for all events was 65°. The average dip was 80°. These 
average values were close the values obtained by fitting of the fault plane by 
least squares and composite focal mechanisms (FPFIT). 

To estimate SHmax with individual focal mechanism, we did an average of 
the P-axes directions (for the events of mechanism strike-slip) and B-axes (for 
the events of normal mechanism). The obtained value was 258°, between the 
values obtained by FPFIT for direction of the P axis (278°) and the value of 
SHmax for the northwestest Ceara (293°) obtained by Ferreira et al. (1998). 

6 Discussion 

A network composed by seven portable digital stations was used to mon- 
itoring the seismic activity in the Tucunduba Dam, almost all installed in 
granitic/gneissic bedrocks of the crystalline basement. It was possible to ob- 
tain records of good quality, with clear P and S waves. Besides, the network 
configuration has two stations that were installed close to the epicentral area. 
Thus, it was possible to determine the hypocenters with good accuracy. 

The active zone was too small, less than about 1 km length and depth 
varying from 4.5 to 5.2 km, considering all events analyzed. It was inside the 
Tucunduba Reservoir. 
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The set of 16 events recorded by the same six stations was selected from 
the method of the least square, obtaing values of 60° for the azimuth and 88° 
for the dip. 

The focal mechanism solution shows a NE-SW fault, with strike-slip dex- 
tral movement with normal component. The fault plane was chosen by taking 
into account the hypocentral distribution (Figure 6). The direction of the 
axes P (278°), obtained by composite focal mechanism, differs 15° of the 
direction of the horizontal maximum stress (Snmax), obtained by Ferreira et 
al. (1998). In the determination of individual focal mechanisms, there is 
three or four stations with the same polarities, which may be indicated that 
mechanisms is not change, thus the strike-slip mechanism (four) is fully with 
expected results. The average direction of the individual planes ranges from 
65° to the azimuth and 80° for the dip, which are close to the values obtained 
by the composite mechanism and least squares. As the direction of horizon- 
tal maximum stress (Snmax) was 258°, it differs 6° of the best estimated by 
Ferreira et al. (1998). 

The seismic activity is located in the border of the basin, common charac- 
teristic of the seismicity of the Northeast (Assumpgao, 1998; Ferreira et al, 
1998). The focal mechanism and fault plane estimated (the strike-slip type) 
is predominant in the Northeast of Brazil (Assumpgao et al., 1985; 1989; 
Ferreira et al, 1995; 1998). The Sumax estimated has a small difference. 
Probably, that difference is because Ferreira et al, (1998) used events to the 
South of the active area and that the line of the coast suffers more flexural 
effect than area to the South. 

The NW area of Ceara was analyzed up to now with events that happened 
on the basements, the seismic area of the Tucunduba reservoir is consistent 
with stress regime (Ferreira et al, 1998), in other words, Sumax parallel the 
Northeast coast, where this active is of the border of the boundary basin. 
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Table 1: Events recorded by the same the six stations (SN01, SN02, SN04, 
SN05, SN06, SN07 and SN09). H-origin is the origin hour, mo is the dura- 
tion magnitude using the same parameters presented by Blum & Assumpgao 
(1990) (m D = l.QAlogD - 0.97). 



DATA H-origin lat S long W depth mo 



1 


970904 


2052 


47. 


.87 


3- 


•11. 


.63 


40-25. 


.79 


4, 


.87 


1.5 


2 


970905 


0113 


50. 


93 


3- 


•11. 


,64 


40-25. 


.81 


4. 


.81 


1.0 


3 


970905 


0114 


02. 


.56 


3- 


•11. 


,63 


40-25. 


.78 


4. 


.80 


1.0 


4 


970909 


0427 


55. 


.67 


3- 


•11 


.70 


40-25. 


.88 


4. 


.78 


1.5 


5 


970911 


2146 


59. 


.60 


3- 


■11 


.66 


40-25. 


.88 


4. 


.70 


< 1.0 


6 


970912 


0517 


31 


.43 


3- 


•11. 


.66 


40-25. 


.87 


4. 


.74 


< 1.0 


7 


970912 


0517 


59. 


.47 


3- 


■11 


.64 


40-25. 


.84 


4. 


.79 


< 1.0 


8 


970912 


0518 


09. 


73 


3- 


■11 


.66 


40-25. 


.86 


4. 


.70 


< 1.0 


9 


970912 


0519 


30 


.63 


3- 


•11. 


.68 


40-25. 


.88 


4. 


.70 


< 1.0 


10 


970912 


0519 


40. 


.12 


3- 


•11. 


.67 


40-25. 


.81 


4, 


.78 


< 1.0 


11 


970912 


0519 


55. 


12 


3- 


•11. 


.69 


40-25. 


.86 


4. 


71 


< 1.0 


12 


970912 


0523 


55. 


.15 


3- 


■11 


,68 


40-25. 


.87 


4. 


.68 


2.2 


13 


970912 


0536 


47. 


.90 


3- 


■11 


.64 


40-25. 


.81 


4. 


.76 


1.5 


14 


970912 


0548 


26 


.82 


3- 


•11 


.65 


40-25. 


.80 


4. 


.74 


1.5 


15 


970912 


0551 


04. 


.71 


3- 


■11 


,65 


40-25. 


.84 


4. 


.79 


2.3 


16 


970912 


0843 


28. 


.86 


3- 


11 


,68 


40-25. 


.85 


4. 


73 


< 1.0 



Table 2: Result of the composite focal mechanism using FPFIT program 
(Reasemberg & Oppenheimer, 1985). 



Azimuth dip 


Rake 


P-Az/plunge 


T-Az/plunge 


60° ±11° 65° ±5° 


-174° ±4° 


278° 27° 


15° 14° 
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Table 3: Results of focal mechanism - FOCMEC. Each solutions with its 
mechanism focal and teh events correspond Tabled 



Events 


Azimuth 


dip 


Rake 


Mechanism 


Az. /plunge 
eixo-P 


Az. /plunge 
eixo-T 


1 


44° 


68° 


-166° 


strike-slip 


265°/25° 


358°/06° 


6 


74° 


44° 


-84° 


normal 


84°/86° 


340°/01° 


11 


78° 


85° 


174° 


strike-slip 


125°/00° 


35°/08° 


13 


76° 


65° 


-45° 


normal + strike-slip 


215°/49° 


315°/8° 


14 


89° 


84° 


-10° 


strike-slip 


45°/12° 


135°/03° 


15 


51° 


72° 


-172° 


strike-slip 


263°/33° 


360 o /H° 


16 


44° 


73° 


170° 


strike-slip 


89°/07° 


358°/07° 
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Figure 1: Map of the locations of the seismographic stations. The triangles rep- 
resent the seismographic stations. In the site SN01 were deployed an analogical 
(SN1A) and a digital station. The squares represent the small villages; lines cross- 
ing the dam are faults (DNPM, 1973). 
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Figure 2: Map separating the two main domains NW of Ceara. The Medium 
Coreavi Domain = MCO, the Ceara Central Domain = CC and transbrasiliano 
lineament. Cicles are epicenters recorded by UFRN network. 



15 



200 



150 

o 

03 100 



05 
Q 



50 







J i I i I i I i I i I i I i I i I i L 



-Jrm b = 3.0 



40 80 120 160 200 240 280 320 360 400 

Jun/1 1/1997 to Aug/14/1998 



Figure 3: Tucunduba activity, Number of daily events recorded by SN1A station 
from June 1997 to August 1998. 
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Figure 4: Seismogram recorded at SN07 station with clear P and S arrivals in the 
vertical component. Event n° 15 (table 
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Figure 5: Map of Tucunduba dam epicenters, with 160 selected events (circles), 
recorded by at least four stations (rms < 0.02s, erh < 0.02 km, erz < 0.2 km). 
Triangles indicate seismic stations. 
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Figure 6: Map of epicenters with 16 selected events (circles), recorded by at least 
six stations. Triangles indicate seismic station. Small squares represented a closer 
look of the picenters. 
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Figure 7: A. Detailed map of epicenters in Figure 6. B e C. Projections of the 
hypocenters on the vertical plane across the probable fault plane (Parallel and 
Perpendicular), respectively. Events used to estimate fault plane. 
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Figure 9: Composite Focal Mechanism for 16 events (table 1). Lower hemisphere, 
equal-area projections. Black and Ash Circles denote compressive and extensive 
first motions dominated, respectively. 
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